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Fast evolving size of early-type galaxies at z > 2 and the 
role of dissipationless (dry) merging 
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ABSTRACT 

Wc present the analysis of a large sample of early-type galaxies (ETGs) at < z < 3 
aimed at tracing the cosmic evolution of their size and compare it with a model of 
pure dissipationless (dry) merging in the ACDM framework. The effective radius i?e 
depends on stellar mass Ai as Re{A4) oc with a ^ 0.5 at all redshifts. The redshift 
evolution of the mass- or SDSS-normalized size can be reproduced as oc (1 -I- z)^ with 
/3 ~ — 1, with the most massive ETGs possibly showing the fastest evolutionary rate 
(/3 —1.4). This size evolution slows down signifieantly to (3 ^ —0.6 if the ETGs at 
z > 2 are removed from the sample, suggesting an accelerated increase of the typical 
sizes at z > 2, especially for the ETGs with the largest masses. A pure dry merging 
ACDM model is marginally consistent with the average size evolution at < z < 1.7, 
but predicts descendants too compact for z > 2 progenitor ETGs. This opens the 
crucial question on what physical mechanism can explain the accelerated evolution at 
z > 2, or whether an unclear observational bias is partly responsible for that. 
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1 INTRODUCTION 

Early-type galaxies (ETGs) are important probes of struc- 
ture formation and massive galaxy evolution. At < z < 1, 
the ETG stellar mass function shows a downsizing evolution 
apparently difficult to reproduce with the current models 
of galaxy formation, with the majority of massive ETGs 
{M > 10^^ Mq) already in place at z « 0.7, (Pozzetti et 
al. 2010 and references therein). At z > 1, the information 
is still incomplete, but bona fide ETGs have been identified 
up to z ~ 2.5 (e.g. Kriek et al. 2006; Cimatti et al. 2008 
and references therein). These high-z ETGs are character- 
ized by old stars (1-3 Gyr), e-folding decaying star forma- 
tion timescales r ~0. 1-0.3 Gyr, low specific star formation 
rates (SSFR< lO^'^ Gyr^^), low dust extinction, large stellar 
masses (A^ > 10^ '^Mq), spheroidal morphologies (although 
some of these systems have a disk- like component; van der 
Wei et al. 2011), and number densities growing rapidly from 
z > 3 to z ~ 1 (e.g. Fontana et al. 2009; Dommguez-Sanchez 
et al. 2011; Brammer et al. 2011). 

A puzzling property of ETGs at z > 1 is that they have 
smaller sizes, down to effective radii i?e < 1 kpc, and corre- 
spondingly higher internal mass densities than present-day 



ETGs with the same mass (e.g. Daddi et al. 2005; Trujillo 
et al. 2006; van der Wei et al. 2008; Cimatti et al. 2008; 
Buitrago et al. 2008; Saracco, Longhetti & Andreon 2009; 
Williams et al. 2010; Cassata et al. 2011; Damjanov et al. 
2011; Newman et al. 2011, and references therein). It is not 
clear yet whether the environment plays (e.g. Papovich et al. 
2011; Cooper et al. 2011) or not (e.g. Rettura et al. 2010) 
a role in the ETG size evolution. The few available mea- 
surements of stellar velocity dispersions are consistent with 
those expected from the ETG sizes and confirm that these 
systems are truly massive (Cenarro et al. 2009; Cappellari et 
al. 2009; Onodera et al. 2010; van Dokkum, Kriek & Franx 
2009; van de Sande et al. 2011). Several models have been 
proposed to explain the size-mass evolution, including dissi- 
pationless (dry) major and minor merging, adiabatic expan- 
sion driven by stellar mass loss and/or strong feedback, and 
smooth stellar accretion (e.g. Khochfar & Silk 2006; Bour- 
naud. Jog & Combes 2007; Fan et al. 2008; Nipoti et al. 
2003, 2009a, 2009b; Naab et al. 2009; Hopkins et al. 2009; 
Oser et al. 2011). However, the global picture is far from be- 
ing clear. In this paper, we exploit a large sample of ETGs in 
order to investigate the evolution of their size as a function 
of redshift and mass, and compare it with the predictions of 
cosmological models of structure formation. 
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Table 1. ETG subsamples 



Sample 


N 


Redshift 


Age 


Ref. 


SDSS 


59500 


< z < 0.4 


yes 


1 


COSMOS/zCOSMOS 


950 


< z < 1 


no 


2 


GOODS-N+S 


469 


< z* < 2 


yes 


3 


Literature 


465 


0.2 < z < 2.7 


no 


4 


GMASS 


45 


1.4 < z* < 3 


no 


5 


COSMOS 


12 


1.4 < z* < 1.8 


yes 


6 


XMMU J2235-2557 


11 


z=1.39 


yes 


7 


K20-0055 


9 


0.7 < z < 1.2 


yes 


8 


POWIR 


6 


1.2 < z* < 1.8 


no 


9 


K20 


4 


1.6 < z < 1.9 


yes 


10 


1255-0 


1 


z = 2.186 


yes 


11 


FW-4871 


1 


z = 1.902 


yes 


12 



z*: a fraction of redshifts is photometric, Age: avail- 
able stellar ages. 1. Hyde & Bernardi 2009, 2. 
|http: //cosmos. astro.caltech.edu/data/index.html[ Scarlata 
et al. 2007, Moresco et al. 2010, 3. Cassata et al. 2011, 4. 
Damjanov et al. 2011, 5. Cassata al. 2008, 6. Mancini et al. 2010, 
7. StrazzuUo et al. 2010, 8. di Serego Alighieri et al. 2005, 9. 
Carrasco, Conselice & Trujillo 2010 (nsersic > 2), 10. Cimatti 
et al. 2004, 11. van Dokkum et al. 2009, 12. van Dokkum & 
Brammer 2010. 

2 THE SAMPLE 

In order to improve statistically on previous studies, we se- 
lected a large sample of 1975 ETGs at 0.2 < z < 3 by 
collecting data from the literature and public data, requir- 
ing the availability of spectroscopic redshifts (or high-quality 
photometric redshifts for z > 1.4), stellar masses, sizes {Re) 
and, when possible, age of the stellar population (see Tab. 
1). The ETGs in the different subsamples share the global 
property to have been originally selected based on the com- 
bination of colors, spectra (or sometimes also SSFR) typical 
of old/passive galaxies with the confirmation of spheroidal 
(E/SO) morphology, or vice versa. We recall that selec- 
tion criteria for ETGs are strongly correlated, with up to 
^85% color/spectra-selected ETGs being also morphologi- 
cally E/SO (e.g. Renzini 2006 and references therein). The 
ETG sizes were generally measured in the observed-frame 
red-optical for low-/intermediate-redshift samples and/or in 
the near-infrared for higher redshifts, i.e. typically sampling 
the rest-frame optical region at all redshifts. Recently, Dam- 
janov et al. (2011) and Cassata et al. (2011) have shown that 
the sizes measured in the rest-frame UV and in the optical 
correlate very strongly with each other, thus excluding sub- 
stantial biases dependent on the wavelength at which the 
size was measured. The SDSS sample of Hyde & Bernardi 
(2009) was included as the reference sample at 2 ~ 0. 

The different subsamples were harmonized to the same 
cosmology [Ho = 70 km s"^ Mpc~\ Q,^ = 0.3, Q.a = 0.7), 
and the stellar masses and ages rescaled to the Maraston 
(2005) stellar population synthesis models with the Chabrier 
Initial Mass Function (IMF) by using empirical scaling re- 
lations of Pforr et al. (2012). 

Based on the information available in the litera- 
ture, the stellar mass completeness of most subsamples is 
log A^comp/Af© = 10.5 at all redshifts. However, a few sub- 
samples have logXcomp/M© ~ 10.8 - 11: GN/DEIMOS 
(from Damjanov et al. 2011) and zCOSMOS aX z < 1.3, 
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Figure 1. The size - stellar mass relation for log A^/Mq >10.5. 
Black and blue symbols indicate SDSS and non-SDSS ETGs re- 
spectively. Only SDSS ETGs older than 10 Gyr are shown in order 
to avoid excessive crowding. The black solid line shows the size- 
mass relation of Shen et al. (2003) and its ±lcr scatter (dashed 
black lines). Red dashed lines: best fit relations using Re oc Al". 



MUNICS, K20, Mancini et al. (2010) and van Dokkum et al. 
(2008) at 1.3 < z < 2, Cassata et al. (2011) and van Dokkum 
et al. (2008) at z > 2. The final sample was reduced to 1080 
galaxies in order to avoid mass incompleteness effects (see 
Section 4 for details). 



3 THE SIZE - MASS RELATION 

Figure 1 shows the size-mass relation in three redshift 
ranges with a comparable number of galaxy in each bin. A 
power-law fit, Re (x Ai" , applied to non-SDSS ETGs with 
log M /Mq >10.5 provides a = 0.52±0.05, 0.47±0.04, 0.50± 
0.04 from low to high redshift in the three bins of Fig. 1, with 
Q basically independent of redshift. For instance, the ETGs 
with z > 2 have a = 0.45 ± 0.11. This result is consistent 
with recent works (Damjanov et al. 2011; Newman et al. 
2011), and does not depend significantly on the choice of 
the redshift bin limits or the minimum stellar mass cuts. In 
comparison, the SDSS ETGs with logA^/A/© > 10.5 have 
a =0.58±0.01, consistent with Shen et al. (2003). 



4 THE SIZE - REDSHIFT RELATION 

Fig. 2 shows the redshift evolution of the ETG size 
in two complementary ways: the mass-normalized radius 
(i?e(2)/A4fi, where Mii = X/IO" Mq, adopting a = 0.55 
as representative value), and radius normalized to the aver- 
age size of SDSS ETGs {Re{z)/Re[SDSS)) in three mass 
bins. Both quantities are useful to derive the size evolu- 
tion independently of the correlation between Re and M. 
The evolution is parametrized by the usual functional form 
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Figure 2. Bottom panel: the mass-normalized radius (Re/Aljj^'^) 
of ETGs with log M/Mq >10.5 as a function of z. Top three 
panels: the fractional radius (relative to SDSS; {Re{SDSS)) = 
2.51, 3.10, 5.61 kpc). In all panels, green and blue small sym- 
bols indicate individual SDSS and non-SDSS ETGs respectively, 
whereas the large black and red filled circles arc the average val- 
ues {(Re/Re{SDSS))) of SDSS and non-SDSS respectively, each 
with the ±1(T scatter. The black solid line is the best fitting func- 
tion to the red and black points. The redshift bins adopted for 
the non-SDSS ETGs in this figure are < z < 0.5, 0.5 < 2 < 1, 
1 < z < 1.5, 1.5 < 2 < 2, 2 < 2 < 3. The cyan dashed Unes 
show the best fits for the case of ETGs with 2 < 2 (redshift bins: 
< 2 < 0.3, 0.3 < 2 < 0.6, 0.6 < 2 < 1, 1 < 2 < 1.5, 1.5 < 2 < 2). 



size oc (l + z)^ . Clearly, this parametrization does not mean 
that all high- 2: ETGs are the direct progenitors of all low- 
z ETGs because these galaxies evolve in the redshift range 
< z < 3 through several processes and increase their num- 
ber density and mass, but it has simply the statistical mean- 
ing of indicating how the typical sizes compare at different 
redshift s. 

In order to mitigate the potential effects of stellar mass 
incompleteness, for each mass bin of Fig.2 (top three panels), 
the galaxies with M < Aicomp{z) were removed from each 
subsample. Thus, each mass bin is always complete down to 
the minimum mass of the bin. For instance, the ETGs with 
M < 10" Mq at z > 1 and < 10^° '^ Mg at z > 2 have 
been excluded from the zCOSMOS and the Cassata et al. 
(2011) sample respectively. 

For the mass- normalized radius (Fig.2, bottom panel), 
we derive /3= -1.06 ±0.14(-1.24 ±0.15) for logX/M© > 
10.5(10.9), in agreement with the literature (Damjanov et 
al. 2011, Newman et al. 2011). This result is stable against 
changing a between 0.4 and 0.7, the boundaries of the red- 
shift bins, and the minimum stellar mass. The consistency of 
the Re (X M" and (Re/Mii) oc (1 -|- z)^ relations with pre- 
vious results confirms that our sample, despite its somehow 
heterogeneous composition, can be reliably used to perform 
further detailed studies. 

Fig. 2 (top three panels) shows Re normalized to 



Re{SDSS). For the three bins of increasing mass, we de- 
rive /3=-1.23± 0.15, -0.87±0.14 and -1.39±0.13. This sug- 
gests that /3 does not significantly depend on stellar mass, 
although a steepening is suggested in the most massive bin. 
These results are stable, within the statistical uncertainties, 
against changes in the minimum stellar mass in the lowest 
mass bin, in the boundaries of the mass bins, and in the 
minimum age of the SDSS ETGs. 

Fig. 2 also shows that the ETG size growth rate seems 
to increase significantly at high redshift, suggesting a slower 
evolution at z < 2. For the mass-normalized radius (bot- 
tom panel), /3 becomes -0.57±0.15 if the ETGs at Zcut > 2 
are excluded. Similarly for Re/Re{SDSS), P becomes /3=- 
0.73±0.14, -0.76±0.13 and -0.75± 0.14 if ETGs &t z > 2 are 
excluded in the three bins of increasing mass. We note that 
the flattening is particularly pronounced and significant in 
the bin of largest masses. The flattening of /3 is also present, 
within the statistical uncertainties, against changing the lim- 
its of the redshift bins if we cut ETGs around 1.5 < Zcut < 2. 
For instance, if the ETGs with Zeut > 1.5(1.7) are ex- 
cluded in the bin with log M/Mq >10.9, we derive /3=- 
0.95(-0.72)±0.15. The overall results presented in this sec- 
tion do not change if median values are used instead of the 
average ones, nor if the SDSS data points are excluded. 

The potential role of the so-called progenitor bias (e.g. 
Saglia et al. 2010 and references therein) has been assessed 
through an age filtering by comparing the size of ETGs hav- 
ing ages compatible with the cosmic time passed from high- 2: 
to low-z. For a given redshift range z\ < z < Z2 and aver- 
age redshift z, we estimated /3 by comparing the average 
size of the ETGs (having an average age tage(z) ± <Jt(z)) 
with the size of SDSS ETGs with ages at zq between 

[tage{z) — Cri(.)] -f t{z - Zq) aud \tage{z) + 0-i(-)] + t{z - Zq) , 

where t[z — zo) is the cosmic time passed from z to zo. Based 
on this approach, no significant or systematic changes of /? 
have been found for a variety of tested redshift ranges. For 
example, for ETGs with log A^/A/q >10.9 at 0.7 < z < 1.3, 
we derive /? = -1.06 ± 0.2 and /3 = -0.97 ± 0.2 respectively 
with and without applying the above age filtering. Similarly, 
if we select the most distant ETGs (1.8 < 2 < 3), we obtain 
/3 = -1.29±0.2 and/3 = -1.23±0.2. If we take these results 
at face value, this implies that most of the Re{z) evolution 
is unlikely to be the result of a progenitor bias due to high-jz 
ETGs being preferentially selected to be redder and more 
compact than lower-2 younger and larger ETGs missed at 
high-2. However, we recall that, due to the heterogeneous es- 
timates of the stellar ages in our sample, this result should 
be confirmed with larger and more homogeneous samples. 
On the other hand, we note that large ETGs are indeed 
absent at 1.4 < 2 < 3 in the GMASS subsample (Fig. 2), 
that is selected based solely on morphology irrespective of 
colours (Cassata et al. 2008). 



5 COMPARISON WITH A ACDM MODEL 

Dissipationless (dry) merging is one of the few mechanisms 
known to make galaxies less compact (e.g., Nipoti et al. 2003; 
Naab et al. 2009), so it is often invoked to explain the ob- 
served size evolution of ETGs. Here we briefly address the 
question of whether the observed size evolution is consistent 
with the merger histories of concordance A cold dark matter 
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(ACDM) cosmology. For this purpose, we use the ACDM- 
based merger models presented in Nipoti et al. (2012, here- 
after N12), which are such that the variation in surface- 
mass density is maximized, because dissipative effects are 
neglected. In particular, we adopt here model B of N12, 
which is the most strongly evolving of their models, so the 
predicted size evolution must be considered an upper limit. 
Here we briefly describe the main properties of the model, 
but we refer the reader to N12 for details. For an observed 
ETG with measured stellar mass M and effective radius 
Re, the model allows to calculate the redshift evolution of 
M and Re via analytic functions calibrated on A'^-body sim- 
ulations. The galaxy is first assigned a halo mass using the 
redshift-dependent stellar to halo mass relation of Behroozi, 
Conroy, & Wechsler (2010). The halo growth- history is then 
computed using Fakhouri et al. (2010) fit to halo merger 
histories in the Millenium I and II simulations (Springel et 
al. 2005; Boylan-Kolchin et al. 2009). The associated growth 
of M is obtained by assigning stellar mass to satellite ha- 
los with Behroozi et al. (2010) recipe, and considering only 
mergers with mass ratio >0.03 (to exclude cases with too 
long merging time). Finally, the corresponding variation in 
Re is computed using analytic functions verified with A^- 
body simulations of minor and major dry mergers between 
spheroids. For observed high- 2: ETGs, we compute the pre- 
dicted M{z) and Re{z) up to z = 0.14, which is the average 
redshift of massive (logM/MQ > 10.9) SDSS ETGs. The 
evolution in the M-Re plane is shown in Fig. 3, taking as 
progenitors the observed ETGs with 2 < 2: < 3 ((2:) ~ 2.4) 
and those with 1.5 < 2 < 2 ((2) ~ 1.7). In both cases the 
present-day descendants tend to be more compact than real 
2 ~ ETGs, but the deviation from the SDSS M-Re rela- 
tion is larger than the observed scatter ( ~ 0.15 logi?e at 
given M) only when 2 ~ 2.4 progenitors are considered: the 
2 ~ descendants have median vertical offset from the SDSS 
best-fit Alogiie ^ -0.3(-0.1) for 2 ~ 2.4(1.7) progenitors. 
In the case of 2 ~ 2.4 progenitors, also the 2 = 1 model 
descendants are more compact than real 2 ~ 1 ETGs. We 
conclude that the 2 > 2 ETGs are so compact that, even 
according to extreme pure dry-merger models, their low-2 
descendants are predicted to be significantly more compact 
than present-day ETGs. On the other hand, the milder size 
evolution observed since 2 ~ 1.7 is marginally consistent 
with ACDM dry-mergei models, though the model descen- 
dants are distributed in the A4-Re plane with larger scatter 
than the observed ETGs (see also N12). An additional prob- 
lem is the presence, among the predicted 2 ~ descendants, 
of outliers, i.e. galaxies in regions of the M-Re in which there 
are no SDSS ETGs: for instance, three model galaxies with 
log A^/M0 ~ 12 and Re > 70 kpc, and a model galaxy with 
log AI/Mq ~ 11.5 and Re ~ 0.3 kpc (Fig. 3). We recall here 
that the existence of compact low-2 ETGs with sizes and 
masses comparable to those of compact ETGs at 2 > 2 is 
somehow unclear, with some results showing an absence of 
such galaxies (e.g. Taylor et al. 2010) and others finding a 
few candidates (e.g. Shih & Stockton 2011; Valentinuzzi et 
al. 2010). Do our results necessarily imply that dry merging 
alone cannot explain the observed size evolution ? In princi- 
ple it could be the case that dry mergers are responsible for 
the whole size evolution, but the actual rate of mergers is 
higher than predicted. This hypothesis can be tested by fur- 
ther comparison with observations. A first constraint comes 



from the observed redshift evolution of the ETG stellar mass 
function. Let us take, for instance, the model describing the 
evolution of the 2 ~ 1.7 ETGs. At each 2 < 1.7 we se- 
lect only model galaxies with log M/Mq > 10.9 and for this 
subsample we measure the average stellar mass (log M). The 
redshift variation of (log A^) is found to be well represented 
by the fit (logX/Mo)(2) = 11.57 - 0.152, i.e. the aver- 
age mass increase by ~ 30%(70%) from 2 = 0.7(1.5) to the 
present. This is compatible with observed evolution of the 
ETG stellar mass function at < 2 < 1 (e.g. Pozzetti et al. 
2010). Another testable feature of the model is the predicted 
merger rate. Let us define major(minor) mergers those with 
mass ratio>(<)l/4. The predicted number of major merg- 
ers per unit time dN,n/dt decreases for decreasing 2 (see 
Fakhouri et al. 2010): for our model ETGs we get, on av- 
erage, dNm/dt ~ 0.13(0.23) Gyr"^ at 2 ~ 0.55(1.15). These 
rates are higher by a factor of ~ 2 than estimated obser- 
vationally at similar 2 by Bundy et al. (2009; see also Lotz 
et al. 2011), indicating that the considered model might be 
extreme also in this respect. In the model, both major and 
minor mergers contribute significantly to the growth of stel- 
lar mass (for instance, ~ 50% each between 2 = 1.5 and 
2 = 0), so, at least within ACDM, massive ETGs do not ac- 
crete most of their mass in very minor mergers, which would 
be more effective in increasing the galaxy size. The above 
arguments suggest that other processes not included in the 
model should contribute significantly to the size evolution. 
Unfortunately, at the moment the proposals for additional 
mechanisms are not very promising: Fan et al. (2008) envis- 
aged that feedback from QSOs could play a role, but also 
this scenario is not without problems (Ragone-Figueroa & 
Granato 2011). 



6 CONCLUSIONS 

The analysis of a large sample of ETGs at < 2 < 3 shows 
that their size evolves independently of stellar mass and pos- 
sibily faster at 2 > 2 (especially for ETGs with the largest 
masses). The interpretation of this result is not straight- 
forward as the available information does not allow us to 
assess if this is an observational bias (e.g. large ETGs with 
low surface brightness are missed in high-2 samples), or it is 
an intrinsic change in the evolutionary pattern implying a 
very rapid growth of ETGs from 2 > 2 to lower redshifts. We 
explored the possibility of pure dry merging as the dominant 
growth mechanism within the ACDM framework, and found 
that this scenario is marginally consistent with the average 
size evolution at < 2 < 1.7, but predicts descendants too 
compact for 2 > 2 progenitor ETGs. Further studies and 
larger samples of ETGs at 2 > 1.5, which will be obtained 
with future wide-field surveys (e.g. Euclid, Laureijs et al. 
2011), will shed light on these open questions. 
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Figure 3. Redshift evolution of ETGs in the stellar mass-effective 
radius plane, as predicted by the dry-merger ACDM model de- 
scribed in Section 5. The progenitors (circles) are the ETGs ob- 
served at 2 < 2 < 3 {(z) ~ 2.4, upper panel) and those observed 
at 1.5 < z < 2 ((z) ~ 1.7, lower panel). The predicted descen- 
dants at z = 1 and z = 0.14 are represented, respectively, by the 
triangles and the squares. For comparison we plot SDSS ETGs 
(green dots) and the local SDSS best-fit with its observed scatter 
(black solid and dashed curves; Shon et al. 2003). The red dashed 
lines indicate the best-fits for ETGs observed at 0.8 < z < 1.2. 
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